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Abstract 

The apparent nonlocality of the Coulomb gauge external field prob- 
lem in electrodynamics is illustrated with an example in which non- 
locality is especially striking. A rather intuitive derivation of this 
apparent nonlocal behaviour from a purely local picture is given. In 
the course of the derivation Lorentz-force naturally decomposes into 
a field and source part and a novel classical effect based on this sepa- 
ration is pointed out. 



1 Az introductory example 

Consider a rigid body with a given charge distribution p(x) on it but of 
total charge equal to zero. If it is placed into a homogeneous magnetic field 
B = V x A then in the combined electric field E of the body and the 
external magnetic field a certain amount N of angular momentum will arise. 
The simplest way to demonstrate this is to slowly remove the magnetic field 
and to calculate the torque 

K(t) = J rf 3 x[x,pE'(t)] (1) 

acting on the body due to the electric field E'(i) induced by the changing 
B(t). The total amount of angular momentum imparted to the body is equal 
to the angular momentum N. 
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In Coulomb gauge the vector potential is chosen so as to make its diver- 
gence to vanish and the scalar potential is equal to the Coulomb field of 
the charge density p(x). For a homogeneous field we hav^j] 



A = i[B,x] (V-A = 0). 
The induced electric field is determined by the equation 

E'(t) = --A(t), A(0) = A, A(oo) = 0. 

c 

Then ([T]) can be written as 

K(t) = -| J d 3 x[x, ^A(t)]. 



(2) 



Since N = K the total amount of angular momentum received by the body 
is equal to 

N = + J d 3 x[x,^A]. (3) 
The angular momentum stored in the field can, therefore, be written as 



N = __ / rf 3 xp(x) X) [B)X ] 



Since 



x, [B, x] 



i(3(x-B)x-r 2 B) + jjr 2 B, 



where r 2 = x 2 , the Cartesian components of N are 

N a = —^D a pBp + ^-(pr 2 )B 
be 6c 



(4) 



Here 



D a p = J d 3 xp(x)(3x a xp - r 2 5 af3 ) 



1 The condition V • A = would permit us to add to A a gradient of a harmonic 
function but for the calculation of the gauge invariant quantity (TT]) the choice of the gauge 
is dictated solely by convenience. We note in passing that in this static situation there is 
in fact no difference between Coulomb and Lorentz gauges. 
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is the quadrupole momentum tensor of the body and 



(pr ) = J d x p(x)r 

is its scalar part. 

Simple as it seems our example reveals its rather peculiar feature when 
one tries to recalculate N starting from the equation 



N = / d 6 x 

Anc 



[E, Bl 



(5) 



which defines field angular momentum in terms of field strengths. To this 
end we replace the fields in (j5J) by their respective potentials B = V x A and 
E = — and perform partial integrations within a sphere of arbitrary large 
radius R. Though, neglecting surface terms, we arrive at the expression 
(jlj) again, the contribution of the surface terms turns out to be of finite 
magnitude independent of R (see Appendix): 

Surface Terms = -^~D a pBp. (6) 

Hence we obtain that, contrary to @, within any sphere surroundig the body 
the amount of angular momentum is equal to 

N^ = ~D aP Bp + ^{f^)B a . (7) 

Moreover, this momentum is concentrated within the volume of the body 
since, owing to the independence of (JBJ) of R, in any spherical shell outside 
the body the value of the angular momentum is cancelled to zero. 

The conclusion from this apparent contradiction is that the notion of the 
homogeneous magnetic field of infinite extension is in general an unaccept- 
able idealization. In local problems as e.g. in Zeeman effect, or even in the 
derivation of eq. 01]), the concept of a magnetic field of infinite extension 
is a perfectly suitable idealization but when total field momentum or angu- 
lar momentum are to be calculated the sources of B should be taken into 
consideration explicitely. 

Their simple realization is an infinitely long ideal straight solenoid (an 
idealization too but this time a harmless one) of very large circular cross 
section of radius R s with our body situated on its axis. Then an amount 
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of angular momentum given by ([7j) is concentrated within the volume of the 
body while the remaining part of it, the negative of is found within the 
solenoid above and below the sphere of radius R s (see Appendix). The sum 
total is then given by (jl]) as expected. 



2 The problem of apparent nonlocality 

Assume now that the body we are considering is capable to rotate around a 
fixed point of it. Then its Lagrangian is 

1 If 
L = -I^ttMp + - / d 3 xp(x)A(x) • V(x), 

in which I a p and f2 are the inertia tensor and the angular velocity and V = 
[ft, x]. Using ([3]) and the relation A[f2, x] = f2[x, A], L can be written as 

l = h aP n a n p + n • n. 

In the Hamiltonian H = J ■ Q — L 

8L 

Ja — -QQ - — lafiQp + N a . 

From this 

n a = (r 1 ) afi (Jp-Np). 

Substituting this into H we obtain after some rearragements 

H = \{I- l ) a p{J a - N a ).(Jp - Np). (8) 

The canonical variables in this Hamiltonian are the Euler angles which de- 
termine the orientation of our rotator and their conjugate momenta. We do 
not recall the explicit dependence of J and Q a p on these variables. According 
to (Hj) N depends on the Euler angles through the dipole momentum tenzor 
D a p- 

The Hamiltonian is the energy of the rotator. The Lorentz-force by which 
B acts on the rotator does not perform any work on it and so H is equal to its 
kinetic energy. Hence, the difference (J — N) must be equal to the angular 
momentum R of its rotation. The canonical angular momentum J can, 
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therefore, be interpreted as the sum of the rotational angular momentum 
and of the angular momentum stored in the combined electric field of the 
body and external magnetic field. 

In the previous section we have found that, depending on D a p (i.e. on the 
orientation of the body), a well defined fraction of the field angular momen- 
tum is found at a distance arbitrarily far (say light years away) from the body 
at rest. Remember now that the most characteristic feature of the Coulomb 
gauge is that the potential of a time dependent charge distribution p(x, t) 
is equal to its instantaneous Coulomb potential [1]. Hence, when our body is 
allowed to rotate then the potential determining field angular momentum 
rotates together with it. Hence, in spite of the huge distance, the amount of 
this quantity at the far ends of the solenoid varies synchronously with the 
rotation of the body without any delay. Though from a purely computational 
point of view this apparent nonlocal behaviour presents no difficulty a simple 
explanation of it in terms of local interactions alone would certainly improve 
our understanding of this important aspect of electrodynamic^]. For this 
reason in the next section the results arrived at above will be rederived by 
means of explicitely local considerations. 

3 The underlying physical picture 

Consider now a point charge moving arbitrarily along the trajectory x = r(t) 
in a stationary magnetic field B(x) of finite extension. The Hamiltonian is 



Since the Lorentz-force acting on the charge does not perform any work on 
it this H is numerically equal to the kinetic energy mV 2 /2. Hence p = 
mV + (e/c)A(r) where V = r(t) is the instantaneous velocity of the charge. 

Any gauge transformation of A is accompanied by a canonical transfor- 
mation of p so as to make the difference p — (e/c) A unchanged. As it is well 

2 Another problem related to the nonlocal behaviour of the Coulomb gauge potential 
is the fmiteness of signal (light) velocity. A substantial amount of literature has been 
devoted to this subject since W. Heisenberg proposed it to his disciple S. Kikuchi [2] and 
E. Fermi published his study of quantum electrodynamics [3]. A recent discussion of this 
problem in a variety of gauges is found in [4] where references to earlier literature are also 
found. 
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known (see e.g. [5]) in the special case of the Coulomb-gauge the quantity 
(e/c)A is equal to the momentum G stored in the combined electric field of 
the charge and the external field B when the charge is at rest at the point r(t) 
of its trajectory. The field momentum, therefore, seems as rigidly attached 
to the particle as was the case with the field angular momentum N and the 
rotator in the example of the previous sections. At any moment of time t 
and point P in space the electric field contributing to G is determined by 
the location of the charge r(t) at the same moment of time however far away 
from P. Since composite charged bodies are built up of pointlike charged 
constituents the nonlocality found in the case of the rotating body has been 
inherited from this description of the point charge. 

The sole virtue of the example with the rotator in the previous section was 
that, owing to the splitting of N into fractions located within two disjoint 
domains of space, its nonlocality appeared more spectacular than that of 
a point charge without structure. But in order to elucidate the physical 
picture behind this seemingly nonlocal Hamiltonian description it is sufficient 
to confine ourselves to the case of a single point charge. 

The electromagnetic part of the energy-momentum tensor decomposes 
into the sum of three terms according to the power of B in them. Each 
term satisfies a balance equation (continuity equation with sources) of its 
own. This may be verified either by direct computation or by rescaling 
the external field B with a scalar factor k. For an arbitrary value of k the 
electromagnetic energy-momentum tensor becomes equal to the sum k°T^ + 

^T^p + k 2 Tap- Since the sum satisfies a balance equation for any value of k 
the three parts must satify it separately. 

The term of second order, being constant in time, is obviously of no in- 
terest. The term independent of the external field describes radiation and 
radiation reaction (self-interaction) of the charge and obeys a balance equa- 
tion with a source equal to the negative of the self-force. This term is of 
paramount importance of its own right but for the external field problems 
only the mixed part T^j of the tensor is of significance. The sources of the 
balance equation satisfied by this term in an external magnetic field are the 
negative of the Lorentz-force acting on the charge in the external field and 
the minus Lorentz-force density experienced by the sources of this field from 
the side of the moving charge. From our point of view a cardinal property of 
this part of the energy-momentum tensor is the absence of radiation. This 
follows from the finite extension of the external field in space. 
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The force experienced by the charge at moment to does not depend on its 
subsequent motion. Therefore, in calculating this force we can rely on the 
truncated (at t = to) trajectory which is obtained from the true trajectory by 
bringing the charge to a stop at the position ro = r(to)- At this moment the 
true spatial distribution of the momentum density is certainly different from 
the momentum density (l/47rc)[E^, B] of a charge in state of rest at r but 
the total field momentum will tend to the integral 

G (°) = -L f d 3 x[E^\B] (9) 



Atxc „ 

as t — > oo E). The reason is that due to the absence of radiation the mixed 
field momentum emitted by the moving charge before to will be completely 
absorbed in later times by the sources of the external field. 

Let us compare momentum exchange between the particle and the field in 
the two motions truncated at moments of time to and to + alto- The field mo- 
mentum emitted by the particle in the infinitesimal interval (to < t < to+dto) 
consists of two parts: that what remains stored in the form of field momen- 
tum forever and an additional part which will be subsequently absorbed by 
the sources of the external field. Accordingly, the force experienced by the 
particle consists of two contributions: 

F (o) = _^__ F /( 0)) (1Q) 
dto 

where F'(°) is the Lorentz-force experienced by the sources of B when some 
portion of the mixed field momentum emitted at to by the charge is absorbed. 
Since G^ depends on time through the motion of the charge we have 



dGa ^ 



Vp(to)d p GS> t (11) 



dt 

where d denotes derivation with respect to r. 

The force F'^ is related to the momentum difference received by the 
sources of B in the motions truncated at the moments t + dt and t - Hence, 
in order to evaluate this force we must first calculate the momentum G'^ 
received by the sources via the retarded field of the charge. It is equal to the 
time integral of the Lorentz-force experienced by them: 

G'W = ^-jdt-d 3 x [J(x), B'W(x,t)], (12) 



3 The upper index zero on a quantity indicates that it refers to the truncated motion. 
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where J(x) is the current density supporting the external magnetic field. The 
retarded magnetic field of the charged particle is 

B'(°>(x,t) = V x A'< >(x,i). 

The vector potential here is the Lienard-Wiechert potential of the particle's 
current density j ^ : 

A'(°)(x,t) = - [ df ■ d 3 x' ^ ~\ + ^~ Xl/C) j (0) (x', t'). (13) 

C J |X X | 

Since at t — to the particle comes to a stop we can write 

j(°)(x',t')=j(x / ,t , )0(io-t'), 
where 0(i) is the step-function. Then 

r dtAMfrt) = - [ r dt'i(x',t f ). 

J — oo C J |X X | J—oo 

Now, using the relation 

V x (/U) = [V/, U] + /V x U (14) 
we calculate the curl of this function: 

/oo poo 
dtB'(°>(x,t) = Vx / dtA'(°)(x,t) = 
-oo «/ — oo 



C 











vix-x'r 1 , 


/ cft'j(x',t') 






J — oo 



To obtain F'(°) we have to compute the derivative of G'^ as given in ( TTiZj) 
with respect to to. At first we write 

A r d tB'^(x,t) = - [ d 3 x> [Vlx-x'l^jKto)]. 

Gtfo J-oo C J 

The current density of a pointlike charged particle is 

j(x',t ) = eV(t )5(x'-r(t )). 
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Then 

d l " x 
dto 

where 



/DO p 1 

dfB'<°>(x,t) = - [V|x-r(f )|-\ V(to)] = -[V0^,V(to)], (15) 
-co c c 



|x-r(t )| 

is the potential of a charge at x = r(t ). 
Now we can write 

F'(°) = ^ = ^ 2 Jd 3 x [J(x), [V^\ V(t )]] = 
= rf 3 x{(J-V(t o ))V0W-(J-V0W)V(t o )} 

The second integral on the r.h.s. is equal to zero because the stationary 
current J(x) is solenoidal: 



d x J • V 



( / : W 0) V • J = 0. 



In Cartesian coordinates we have 



<9x r 



<9rv 



-<9 Q 



Then 



f ,(o) = -1$, y (j • v(*„)) (o) = -vp{to)d a i y rf 3 x0w^ 

Equation V x B = (47r/c)J permits us to replace J by B: 
1 / d 3 ^0(°)j« = -L / d 3 ^°V x B. 

C 2 47TC ./ 



(16) 



Here we can use (TT4|) again: 



i y d 3 x<t>wjp = ^~ |y [b, V0 ( ° } ] + y av x (^b; 



The second term is transformed into a surface integral and can be dropped. 
In the first term Vo 1 " 1 = — E^ ). Then, using Q), we obtain 

\ [ d 3 x^J = — [d 3 x [EW, B] = G(°'. 
c 2 J Aire J 

Substitution of this expression into (ITS]) leads to the concise form 

F'} 0) = -V p (to)d a Gf. (17) 

From now on we may drop upper indices. Inserting (|T7|) and (TTTT) into ( fTOj) 
we obtain for the force acting on the charge the formula 

F Q = V^Gp - VpdfsG a , (18) 

where V(£) = r(t) and G is equal to the field momentum stored in the 
combined field of a charge resting at r(t). The vector form of this equation 
is 

F = [V, V x G], 

where V denotes now differentiation with respect to r the only coordinate G 
depends on. 

On the other hand, we know that F is the Lorentz-force which acts on 
the charge: 



F = ^[V, B] 



V, V x -A 

c 



We see that (up to the factor e/c) G can be identified with the Coulomb 
gauge vector potential provided the latter can be chosen vanishing together 
with G and this choice is unique (see the example at the end of this section). 
Equation ( ITBl can, therefore, be written also a£| 

F a = ^(y f} d a Af } -Vpd fi A a ). (19) 

At the beginning of the present section we have stressed the apparent 
nonlocality of this gauge. Now we see explicitely that locality is by no means 



4 Since V is independent of r the covariant form of this equation is 

F= - (V(V • A) — (V • V)A) . 



As explained in the text in Coulomb gauge both terms on the r. h. s. have their own 
physical meaning separately. 
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broken: There is in fact no field momentum attached rigidly to the charge 
and in the general case the true mixed field momentum differs from G at any 
moment of time. The flow of the combined field momentum density obeys 
continuity equation in space between the moving charge and the sources of 
the external field and this is what locality means. Moreover, the physical 
picture revealed provides the interpretation of the terms in the formula ( ITS]) 
for the force F experienced by the moving charge. They describe momentum 
exchange of two different types: the first one with the sources of the external 
field and the second one with the field itself. 

An interesting special example is that of a charged particle passing by 
an infinitely long straight ideal solenoid (a classical 'Aharonov-Bohm situa- 
tion'). Since the magnetic field outside the solenoid is zero the two terms in 
(j!8p cancel each other but neither of them vanishes. That means that the 
charge 'catalizes' momentum transfer between the combined field and the 
coil. Outside the solenoid we have 



where S is the cross section of the solenoid. In this domain the two terms in 
f fTSj) are equal to each other and so ffTTl) can be written as F' a = —VpdpG a . 
Since G = (e/c)A, the total momentum imparted to the coil by the moment 
of time t the charge reaches position r is equal to 



When the charge leaves freely to infinity no net momentum exchange occurs. 
But if it is stopped at r a finite amount f[2~Tj) of momentum is transferred 
from the combined field to the coil. I wonder whether this effect may be 
experimentally observed. 

Inside the solenoid the continuation of the external vector potential ( 120]) 
is given by (J2]) and, as it is well known, the field there is homogeneous. 
Inserting this into ( fT9l) we find that, for a point charge revolving in a plan 
perpendicular to solenoid's symmetry axis on a circle around it, the two 
terms on the r. h. s. are equal to each other. Therefore, the reaction to 
the centripetal force acting on the charge is distributed evenly between the 
field and the coil: the force experienced by the coil is the half of the force 
experienced by the revolving point charge. 




(21) 
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Appendix. Derivation of eq. (H) 

Let us denote the volume integral in (J5J) confined to a sphere of radius R 
by N^. Inserting E = — V0 and B = V x A into it we obtain after some 
rearragements 



N (R) 



where e Q( g 7 is the Levi-Civita symbol. 

In the second term the two V symbols have identical indices. By partial 
integration we make them to merge into V 2 : 



-j—taM [ d 3 x {x p V 5 <p ■ V 7 A - V s (x p ■ V s <f) ■ Ay) + V p <p ■ Ay + x p ■ V 2 <p ■ AJ 

47rc JR 

In the first term the index of the V-symbol and that of the vector potential 
are identical. By further partial integration we form V • A out of them: 



^-e a ^ / d 3 x {V s (x p ■ • V 7 A) - V^Ap) + V 7 • A t 
t7rc Jr 



A-kc 

-<f> ■ a^V^V ■ A) - Vs{x p ■ V 5 ■ Ay) + ■ Ay + xp ■ V 2 ■ A,} 

By eq. (j2j) we have V ■ A = 0. Out of the six terms left we have three volume 
terms but two of them cancel for symmetry reasons: 

e Q /3 7 (V 7 • A p + Vrf • Ay) = 0. 

If we insert V 2 = — 47rp into the remaining one we obtain the r.h.s. of flS} 
the integration being confined to the sphere of radius R: 

N<°> = + 

According to (jlj), its value is equal to 

N^ = -±D a ^ + ±-(pr 2 )B 



[ d 3 x 


fx, ^Al 


Jr 


C - 



for any sphere outside the body. 
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The three surface integrals are 

N£> = +^ e «/37 J d * x V 7 (0A^) 

4vrc J R 

All of them must be of the form (r/j/ 'c)D a pBp (i = 1, 2, 3) and our task is to 
calculate the coefficients rji. 
Using Gauss theorem 



L 



(fVv„/= / <is„-f= I ii 2 <m 



<R JR 

where Q denotes solid angle, we have 

1 R f 

= +-e«/?7^ f d ^ x ^ A fi- (A- 1 ) 

The body is uncharged by hypotheses and so no Coulomb potential con- 
tributes to <fi. The contribution of its dipole moment d to is of the form 
constant x [d, B] but the constant must be zero because it is a polar vector 
while is an axial one. From moments higher than quadrupole no vector 
can be formed with B. Hence it is only the quadrupole potential 

= ■^jpX a xpD a p (A.2) 

which contributes to the surface integrals. 
For the vector potential we can substitute 

Ap = —-ep pa x p B a , 

and we obtain 

= -Y c e aM e Ppa B a-^ j d ^ x y x p (f> = j 5 «^ 2 J dVt (j) + j dVL x c 

The first integral vanishes because of the zero net charge condition. Substi- 
tuting (A.2) into the second term and using the relation 

AttR 4 

dVl X^X nX v X a — (SryuSva -\- S^ v 5ua ^■ya^av) 

10 
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we arrive at the value rji = 1/30. 



A completely analogous derivation leads to r] 2 = —1/30. 



The remaining surface integral is 



A^ 3) = +-e Q/37 — / dQ xpx s V s <j) ■ A 




Since the potential (A. 2) is of degree —3 in r we have 



X S V S (j) = r TT~ = -30. 



Hence 




Comparing this with (A.l) we obtain 773 = +1/10. Adding together N^ ) and 
the surface contributions (j7|) is obtained. 

When the body is situated on the axis of an infinitely long straight ideal 
solenoid of radius R s than (JTj) remains valid for the field angular momentum 
within the sphere provided R < R s . The total field angular momentum, 
however, is given by (jl]) (i.e. N^ )) because the surface integrals vanish as 
R — > 00. Outside the solenoid the vector potential is given by (120]) rather 
than 02]) and tends to zero as R increases. Within the solenoid (EJ) remains 
true at arbitrarily large distances from the body but the surface of the parts 
of the sphere remain finite and the integral vanishes on them for this reason. 

References 

[1] J. D. Jackson, Classical Electrodynamics, 3rd ed. (Wiley, New York, 
1999). 

[2] S Kikuchi, Zs.f.Phys. 66, 558, (1930) 

[http: / / www, springer link, com / content/ v9626 7q8 182ru807/ 

[3] E. Fermi, Rev. Mod. Phys. 4, 87, (1932) 

|http://link.aps.org/doi/10.1103/RevModPhys.4.87| 

[4] J. D. Jackson, Am. J. Phys. 70, 917 (2 002) 
http: / / arxiv.org/abs/physics/0204034 

[5] N. Kroll in Quantum Optics and Electronics, Lectures delivered at Les 
Houches 1964 Gordon and Breach 1965 



14 



